Fast nuclear spin hyperpolarization of phosphorus in silicon 
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We experimentally demonstrate a method for obtaining nuclear spin hyperpolarization, that is, 
polarization significantly in excess of that expected for a thermal equilibrium. By exploiting a mod- 
ified Overhauser process, we obtain more than 68% nuclear anti-polarization of phosphorus donors 
in silicon. This polarization is reached with a time constant of ~ 150 seconds, at a temperature of 
1.37 K and a magnetic field of 8.5 T. The ability to obtain such large polarizations is discussed with 
regards to its significance for quantum information processing and magnetic resonance imaging. 
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Phosphorus doped crystalline silicon (Si:P) is a model 
system for investigating spin effects in the solid state and 
at the same time is a point defect with great technolog- 
ical importance. Si:P has been used since the beginning 
of the semiconductor industry in the early 1950's for ap- 
plications ranging from the ubiquitous (thin film transis- 
tors) to the conceptual (single electron transistors). The 
ability to hyperpolarize the spins in this material is im- 
portant for a number of its applications. Utilizing the nu- 
clear spin of phosphorus donors as quantum bits 1 -! 2 - relies 
on the ability to obtain a well characterized initial state 3 -, 
which can be obtained by hyperpolarization. Spin polar- 
ized silicon microparticles may also have applications for 
magnetic resonance imaging techniques^, similar to other 
hyperpolarized systems, such as xenon^. Whilst it is rea- 
sonably simple to obtain large electron spin polarization, 
for example by using moderate magnetic fields at liq- 
uid 4 He temperatures, doing the same with nuclear spins 
is difficult due to their much smaller Zeeman splitting. 
There are a number of schemes used to obtain nuclear 
spin polarization in excess of the thermal polarization. 
Dynamic nuclear polarization using off-resonance radi- 
ation has been studied extensively 4 ^. Complex pulses 
or adiabatic passage effects may be used to manipulate 
spin states, leading to large polarizations^. Electrical 
injection of hot carriers has been used to obtain posi- 
tive polarizations^ however this requires electrical con- 
tact to the sample. Optical excitation with linearly po- 
larized sub-bandgap light has given small (~ 2.5%) po- 
larization of 29 Si nuclei in silicon with a natural isotopic 
abundance 1 ^. 

In this letter, we demonstrate anti-polarization of 
phosphorus donor nuclei in silicon of up to P = —68%. 
The scheme used is simple, fast and does not involve 
resonant manipulation of either the nuclear or electronic 
spin. Instead, the relative populations are modified using 
photoexcited carriers, generated using white light, at low 
temperatures (about 4 He temperature) and in magnetic 



fields (<~ 8.5 T) significantly smaller than those required 
to obtain an equivalent thermal nuclear spin polarization. 

Phosphorus in silicon can be described by the spin 
(S = 1/2) of its donor electron that is coupled to the 
spin (I = 1/2) of the 31 P nucleus. This model provides 
a system with four energy levels, as shown in Fig. [T] for 
the presence of strong magnetic fields when the nuclear 
Zeeman splitting exceeds the nuclear to donor electron 
hyperfine interaction. At Bq rj 8.5 T, the donor electron 
Zeeman splitting is AE C w 240 GHz whereas the nuclear 
Zeeman energy is AE n « 147 MHz and the hyperfine 
interaction A = 117 MHz. Figure QJa) shows the rele- 
vant spin relaxation processes that occur in the 31 P donor 
atom. Ti is the rate coefficient associated with longitu- 
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FIG. 1: a) Sketch of the energy levels of the four spin eigen- 
states of a phosphorus donor atom in silicon in presence of 
very high magnetic fields. The dashed arrows indicate al- 
lowed transitions with their respective rate coefficients. I\ 
is for longitudinal relaxation processes, Ta and Yb for scat- 
tering with conduction electrons and Yx for the Overhauser 
flip-flop process, b) Sketch of the change from a thermally 
polarized spin ensemble to a hyperpolarized spin ensemble 
for T rcs >> T S pi n . Note that the Tx and Ti processes act 
continuously (not sequentially as illustrated here). 
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dinal relaxation of the electron magnetization towards 
thermal equilibrium with the crystal lattice at tempera- 
ture T spin . T x is the rate coefficient associated with the 
Overhauser spin relaxation process (a flip-flop) between 
the electron and nuclear spinsii. The dependence of the 
Overhauser rate on temperature and magnetic field has 
been described by Pines ct al. 12 who derived an expres- 
sion 
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where s is the sound velocity of silicon, p is the mass 
density of silicon, 7 a multiplicative factor in the range 
fO to 100, / the nuclear spin and A the hyperfine con- 
stant while wo = 9^bB is the Larmor frequency of the 
electrons with g and fi b representing the electron Lande- 
factor and Bohr's magneton, respectively and B the ap- 
plied magnetic field. It is important to note that the 
Overhauser relaxation process serves to return the two 
spin populations 712 and 713 to thermal equilibrium with 
the phonon reservoir, with a temperature T rcs , which is 
not necessarily the same as the spin temperature T sp ; n . 
Due to the constant generation of new excess charge car- 
riers by the illumination, a steady state will be estab- 
lished in which a constant density of hot electrons per- 
sists. As these hot electrons cascade towards the lattice 
temperature, they will emit phonons at a constant rate 
and thus T Tes > Tspin- Note that the phonons will also 
increase T S pm, however, this effect is minimal due to the 
thermal mass of the silicon, which is held constant by 
the helium bath. Differences between T ICS and T sp j n have 
previously been demonstrated using electrical injection 
of hot carriers^. Additionally, the photoexcited carriers 
may scatter with the bound donor electron o 14 ' 15 , caus- 
ing spin relaxation. We capture this process in our rate 
picture by introducing Ta(Tb), the rate coefficient for 
scattering between spin up (down) free electrons and spin 
down (up) bound electrons. 

Feher has previously discussed the effect of the phonon 
reservoir temperature on the polarization of phosphorus 
in silicon^. If the two characteristic temperatures of 
our system are equal, T res = T sp ; n , then the thermally 
(hardly) polarized equilibrium population distribution is 
obtained. However, forcing T res > T S pin by photoexcita- 
tion of charge carriers, we change the steady state pop- 
ulation distribution. The Overhauser process will try to 
achieve thermal equilibrium between states n-z and 713 at 
a temperature T res , and the longitudinal relaxation pro- 
cess will force states (n\ and 712) and (713 and 714) to 
thermal equilibrium at temperature T sp ; n . See Fig. 1(b) 
for a sketch outlining this process. The result of this sit- 
uation is that the population of ri\ becomes much larger 
than the population of all other states, resulting in a 
net nuclear antipolarization, since P = (" 1+ " 2 )~(" 3 +" 4 ) 
Numerical modelling of this process with realistic values 
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FIG. 2: (a) ESR spectra, measured at T=3K with / rcs = 
240 GHz, with (top) and without (bottom) illumination by 
a mercury discharge lamp. The polarization is determined 
by comparing the areas of the two resonances, obtained by 
fitting the data with two Gaussian lineshapes separated by the 
phosphorus hyperfine splitting, AB = 4.17 mT (solid line). 
The small magnetic field offset between the two spectra is due 
to the random offset caused by the superconducting magnet, 
(b) The 31 P nuclear polarization measured as a function of 
time as obtained from EPR spectra, at T — 3K. The solid 
line represents a single exponential decay function which was 
fit to the data. 
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res and T\ indicate that polarization near 
100% is achievable. 

To investigate this effect, we have undertaken elec- 
tron spin resonance (ESR) and electrically detected 
magnetic resonance (EDMR) experiments at B w 8.5 
T, corresponding to a resonant frequency, / = 240 
GH a 16 ' 17 . Similar experiments have been described by 
us elsewher o 18 ' 19 . The samples used in this study were 
similar to those described in reference^. They consist 
of crystalline silicon with (111) surface orientation with 
a phosphorus doping density [P] ~ 10 15 cm~ 3 , with alu- 
minum surface contacts to allow EDMR. 

Figure [2^a) shows two ESR spectra recorded at B k 
8.5 T and T — 3K. The spectra were recorded by sweep- 
ing B through the expected resonance fields. We fit the 
two observed resonances with two gaussian line shapes. 
We can be sure the signal is from phosphorus donor elec- 
trons due to both the g-factor and hyperfine splitting of 
4.17mT. The low-field (high-field) resonance is due to nu- 
clear spins aligned (anti-aligned) with the external field, 
which we will call spin up, | (spin down, j). The res- 
onances are saturated due to the long relaxation times, 
however, we assume that the relaxation times are the 
same and, as a result, can take the area of the resonance 
as a measure of the number of spins that contribute to 
it£. We thus determine the polarization of the sample, 
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P = ■ The lower spectrum was recorded in the dark, 
and shows a nuclear polarization P = —0.008 ± 0.004. 
Next, light from a mercury discharge lamp was shone 
onto the top side of the sample through an optical fibre, 
and the ESR spectra was remeasured (upper spectrum). 
Again, two resonances are visible, however, they have 
different intensities. Here, we determine the nuclear spin 
polarization P = —0.129 ± 0.002. This is a change in 
polarization over the expected thermal polarization by 
a factor 77 = P/Po w — 78. A similar result is obtained 
sweeping Bq in the opposite direction, indicating that the 
polarization is not a passage effect^. 

The polarization model discussed above predicts that 
the time taken to reach a steady-state polarization 
should be limited by the Overhauser rate, since 1/Tx = 
Tx << r^r^r^. By using previously measured 20 low 
magnetic field (B w 340 mT) values for Tx, and extrapo- 
lating to the field used in the experiments presented here 
using Equ. [T] we obtain for the Overhauser time Tx ~ 65 
s, for T rcs = 3K and to = 240 GHz. Figure EJb) shows 
the polarization measured via ESR after light was ap- 
plied to the sample. The data shows a gradual approach 
to a non-equilibrium steady state. The fit of these data 
with a single exponential decay function shows excellent 
agreement and yields a time constant of r = 150 ± 20 s. 
We believe this is in very good agreement with the predic- 
tions of the Overhauser rate made by Pines et al.— , given 
the uncertainty of the low field value (~ 30 hours) at a 
higher donor density, and the extrapolation over nearly 
two orders of magnitude of the magnetic field on which 
the Overhauser rate depends quadratically. 

One aspect of the experiment above suggests that the 
polarization measured with ESR poses a lower limit on 
the maximum polarization obtained. ESR measures the 
polarization in the entire sample; however, only the sur- 
face is illuminated. We expect that, whilst the charge 
carriers will diffuse throughout the sample, they will thcr- 
malize while they diffuse. This will lead to a strong depth 
inhomogeneity of the reservoir temperature and hence a 
depth dependence of the polarization. While the polar- 
ization will be biggest near the surface which is being 
illuminated it will be minimized on the opposite sample 
surface. 

EDMR is a magnetic resonance detection scheme 
which is sensitive to spins close to the illuminated sam- 
ple surface. EDMR relies on the current through a sam- 
ple being influenced by the observed spin state. In Si:P 
at high magnetic fields, we have showni^ that EDMR 
is able to be observed due to to the spin dependent 
capture/emission mechanism described by Thornton and 
Honig 21 , which we have included in our polarization 
model with Ta and Tb- The effect of this process is to 
decrease the current through the sample when resonant 
excitation of the donor electrons occurs. To measure 
EDMR, we thus require free charge carriers, which are 
provided by the illumination used to polarize the nuclear 
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FIG. 3: a) Electrically detected magnetic resonance spectrum 
of Si:P at T = 1.37 K under illumination. The spectrum was 
fit with two gaussian lineshapes, and the nuclear polariza- 
tion determined by comparing the area of the hyperfine split 
resonances. Here, P — —68 ± 1%. b) The polarization deter- 
mined from EDMR measurements for a number of different 
temperatures. The line is a numerical simulation of the ex- 
pected polarization using the rate model presented in this 
paper. The shaded region indicates where our model does 
not hold. See text for details. The measurements in a) and 
b) were obtained using a xenon discharge lamp, c) The polar- 
ization measured simultaneously with ESR and EDMR, as a 
function of the illumination intensity of the mercury discharge 
lamp used to generate the photocarriers, measured at T = 3 
K. The lines are guides to the eye, and are scaled by a factor 
of approximately 2.5 between the ESR and the EDMR. 

spins. Figure [3^) shows an EDMR spectrum recorded at 
T = 1.37 K, the lowest temperature we were able to ac- 
cess. The spectrum was measured with illumination by 
a xenon discharge lamp, and a device current, 2sd = 500 
nA. The microwaves were chopped at a frequency of 908 
Hz, and the change in current was recorded with a lock-in 
amplifier. As with the ESR measurements, the spectrum 
is well fit by two gaussian lineshapes separated by the hy- 
perfine splitting. Again, we use the area of the resonances 
as a measure of the population in each nuclear spin state. 
The polarization measured here is P — —68 ± 1%. This 
corresponds to an enhancement over the equilibrium po- 
larization of r\ ~ 190, and to an effective nuclear spin 
temperature of ~ 5 mK. 

EDMR measurements allow the observation of a 31 P 
subensemble with a significantly more homogeneous 
reservoir temperature than the ESR measurements. We 
therefore use EDMR to test some of the qualitative prop- 
erties of the polarization model described, namely, the 
lattice temperature dependence and the illumination in- 
tensity (and hence reservoir temperature) dependence of 
the observed nuclear polarization. Fig.[3j3) shows the 31 P 
polarization as a function of the lattice temperature. It 
is found to increase monotonically below T ss 3 K. Based 
on the rate model presented in Fig. [H we calculated the 
polarization using the measured lattice temperature and 
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a constant reservoir temperature whose value was chosen 
to fit the experimental data. The simulation results are 
also shown in Fig. [3Jd). The best fit of the simulated val- 
ues to the measured values was achieved for T rcs = 2.7K, 
in agreement with the expectation that hyperpolarization 
vanishes when T sp j n w T Tes . Note that there is significant 
discrepancy between the fit and the data for tempera- 
tures above T spin = 2.5 K. While the calculated data 
predicts no polarization, the measured data shows a clear 
hyperpolarization of P=-6% at 3K. This discrepancy can 
be attributed to our assumption of a constant T res used 
in the calculation. Note that T rcs > T sp j n for these exper- 
iments. Hence, the assumption of a constant T rcs = 2.7 
K becomes unrealistic at T sp j n > 2.7 K. 

In order to further test the polarization model we 
changed the excitation spectrum of the excess charge car- 
riers from the xenon lamp used for the acquisition of the 
data in Fig. [3^a) and (b) to a mercury lamp which has a 
higher spectral temperature. For the latter we measured 
polarization with both EDMR and ESR at a constant 
lattice temperature of T sp j n = 3 K. As shown in Figure 
[3J:) , the EDMR spectra recorded with the mercury lamp 
yield a significantly higher polarization of up to P=-24% 
(instead of 6% at T sp i n = 3 K) , independently of the in- 
tensity over a range of almost one order of magnitude. 
As expected, at low intensities, when the excess charge 
carrier densities drop into a range where the Overhauser 
process is dominated by T S pin> the nuclear polarization 
vanishes and equilibrium appears. The polarization mea- 
sured with ESR was consistently ss 45 % of that mea- 
sured with EDMR, confirming again the inhomogeneity 
of the reservoir temperature throughout the sample. 

Note that while we have demonstrated polarization 
above P — —68%, our model predicts the possibility of 
even higher anti-polarization at lower temperatures and 
higher optical excitation rates. The technical simplic- 
ity of this polarization method suggests that it may be 
beneficial for a variety of technical applications. For in- 
stance, silicon microparticles are biologically inert which 
makes them prime candidates as contrast agents for in 
vivo magnetic resonance imaging. We see no obvious rea- 
son why the polarization technique presented above will 
not provide the same level of polarization in micropar- 
ticles as we have demonstrated in bulk material. Given 
room temperature spin lifetimes > 20 minutes for 31 P 
nuclei in a-Si:H, a disordered material with a bigger de- 
fect density and a larger hyperfine interaction than crys- 
talline silicon, we expect polarization lifetimes of over 
an hour for this material, easily allowing implementation 
of such experiments^. Also, the rapid polarization of 
31 P nuclear spins demonstrated may offer an initializa- 
tion mechanism for 31 P in silicon spin qubits. 

In conclusion, the data presented above demonstrates 
that hyper (anti-) polarization of phosphorous donor nu- 
clear spins in crystalline silicon can be achieved rapidly 
(of the order of a few minutes) by irradiation with above 



silicon bandgap light at low temperatures and high mag- 
netic fields. Polarization in excess of 68% was demon- 
strated, and discussed in terms of a model arising from 
the increased reservoir temperature driven by phonon 
emission during thermalization of photoexcited carriers. 
The qualitative predictions of this model for the polar- 
ization dependence on lattice temperature, illumination 
temperature and intensity have been verified and techni- 
cal applications of this effect have been discussed. 
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